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1.1 Introduction 



■^ In this chapter we will study a particular long wavelength limit of Einstein's 

>— 5 equations with a negative cosmological constant m d + 1 dimensions. In 

00 such a limit we find that Einstein's equations reduce to the equations of fluid 

^^ dynamics (relativistic generalizations of the famous Navier-Stokes equations) 

1—1 in d dimensions. While the motivation for our study lies within the AdS/CFT 

"^ correspondence of string theory, the fluid/gravity correspondence stands on 

o, its own and can be viewed as a map between two classic dynamical systems. 

►>. 1.1.1 Prelude: CFT stress tensor dynamics from gravity 

QQ An important consequence of the AdS/CFT correspondence [Ch.AdSCFT] 

t^ is that the dynamics of the stress (-energy- momentum) tensor in a large class 

: of d-dimensional strongly coupled quantum field theories is governed by the 

L^ dynamics of Einstein's equations with negative cosmological constant in d+1 

^__i dimensions. To begin with, we shall try to provide the reader with some 

^—-^ intuition for this statement and argue that searching for a tractable corner 

^ of this connection leads one naturally to the fluid/gravity correspondence. 

k> In its most familiar example, the AdS/CFT correspondence [Ch.AdSCFT] 

^ asserts that SU{N) M = 4 Super Yang-Mills (SYM) theory is dual to Type 

IIB string theory on AdSs x S . It has long been known that in the 't Hooft 
limit, which involves taking N ^ oo keeping the coupling A fixed, the gauge 
theory becomes effectively classical. However, it was widely believed that 
for any non-trivial gauge theory the resulting classical system would be too 
complicated to be tractable. The remarkable observation of Maldacena in 
1997 was that this field theory intuition is spectacularly wrong. Indeed, not 
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2 The fluid/gravity correspondence 

only is the classical system governing A/" = 4 SYM tractable, it is actually a 
well known theory, viz., classical Type IIB string theory. 

Now, even classically, string theory has complicated dynamics; however 
in the strong gauge coupling (A — )• oo) regime, it reduces to the dynamics 
of Type IIB supergravity (by decoupling the massive string states). More 
interestingly, Type IIB supergravity on AdSs x S^ admits several consistent 
truncations. The simplest and most universal of these is the truncation to 
Einstein's equations with negative cosmological constant, 

E^^ = R^^-lRgf,^ + Ag^^ = , A = - „p7 • (^-l) 

(Note that the AdS curvature radius i?AdS can be scaled away by a change of 
units ; we therefore set -RAdS to unity in the rest of this chapter) . Having thus 
motivated the study of the most beautiful equation of physics, namely Ein- 
stein's equations of general relativity, we now confront the question: What 
does this imply for the field theory? 

Recall that according to the AdS/CFT dictionary there is a one-to-one 
map between single particle states in the classical Hilbert space of string 
theory and single-trace operators in the gauge theory. For instance the bulk 
graviton maps to the stress tensor of the boundary theory. Taking the col- 
lection of such single trace operators as a whole, one can try to formulate 
dynamical equations for their quantum expectation values in the field the- 
ory. While this can be done in principle, the resulting system is non-local in 
terms of the intrinsic field theory variables themselves. 

However, because we can associate the quantum operators (and their ex- 
pectation values) of the gauge theory at strong coupling to the classical fields 
of string theory/supergravity, we know that the set of classical equations we 
are looking for are just the local equations of Type IIB supergravity on 
AdSs X S^. This reduction, whilst retaining lots of interesting physics, still 
turns out to be too complicated from the field theory perspective. For one, 
the space of single trace operators is still infinite dimensional (at infinite N) , 
and relatedly attempting to classify the solution space of Type IIB super- 
gravity is a challenging problem. However, the fact that on the string side 



we can reduce the system to (1.1), implies that there is a decoupled sector 
of stress tensor dynamics in A/" = 4 SYM at large Ajj 

Actually, there is an infinite number of conformal gauge theories which 
have a gravitational dual that truncates consistently at the two-derivative 

^ While this is always true in two dimensional field theories, such a decoupling is not generic in 
higher dimensional field theories (in fact it is not true of A/" = 4 Yang Mills at weak coupling), 
and is in itself a surprising and interesting fact about the A/" = 4 dynamics at strong coupling. 
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level to Einstein's equations with a negative cosmological constant; AA = 4 



SYM theory is just a particularly simple member of this class. Thus (1.1) 
describes the universal decoupled dynamics of the stress tensor for an infinite 
number of different gauge theories. In the first part of this chapter we will 
focus on the study of this universal sector. Later we will generalize to the 
study of bulk equations with more fields, thereby obtaining richer dynamics 
at the expense of universality. 

Given this association between the dynamics of quantum field theory 
stress tensors to the dynamics of gravity in negatively curved backgrounds, 
it is natural to ask - can we do more? Can we for instance classify all possible 
behaviors of stress tensors? On the gravity side we would have to classify all 



possible solutions to ( 1.1 ); this is a laudable goal and various chapters in this 
book are aimed at addressing this question using different approaches. We 
are going to focus on one that naturally follows from the basic organizing 
principle of physics: separation of scales. 

It is well known that in many situations in physics (as well as chemistry, 
biology, etc.), complicated UV dynamics results in relatively simple IR dy- 
namics. Perhaps the first systematic exposition of this ubiquitous fact was in 
the context of finite temperature physics. It has been known for almost 200 
years now that the dynamics of nearly equilibrated systems at high enough 
temperature may be described by an effective theory called hydrodynamics. 
The key dynamical equation of hydrodynamics is the conservation of the 
stress tensor 

Var'^^ = o, (1.2) 

where Vq is the covariant derivative compatible with the background metric 
7(jb on which this fiuid lives. As this equation is an autonomous dynamical 
system involving just the stress tensor, it should lie within the sector of 
universal decoupled stress tensor dynamics. 

Given that the AdS/CFT correspondence asserts that this universal sector 



is governed by (1.1), we are led to conclude that (1.1) must, in an appropri- 
ate high temperature and long distance limit which we refer to as the long 
wavelength regime, reduce to the equations of d-dimensional hydrodynam- 
ics. Indeed, this expectation has been independently verified in [1] and the 
resulting map between gravity and fiuid dynamics has come to be known 
as the fluid/gravity correspondence. In particular, the specific fiuid dynami- 
cal equations, dual to long wavelength gravity in the universal sector, have 



been determined up to the second order in a gradient expansion (cf . ^1.4.4 ) 



Given any solution to the these fiuid dynamical equations, the fiuid/gravity 



map explicitly determines a solution to Einstein's equations (1.1) to the ap- 
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propriate order in the derivative expansion. The solutions in gravity are 
simply inhomogeneous, time-dependent black holes, with slowly varying but 
otherwise generic horizon profiles. 

The main focus of the present chapter is to explain and present the 
fluid/gravity map at the full non-linear level following [1] and subsequent 
work. The connection between these two systems was established and exten- 
sively studied much earlier at the linearized level in the AdS/CFT context 
(following the seminal work [2]). The first hints of the connection between 
fluid dynamics and gravity at the non-linear level were obtained in attempts 
to construct non-linear solutions dual to a particular boost invariant flow 
[3], which provided inspiration for the fluid/gravity map. Such a map was 
also suggested by the observation that the properties of large rotating black 
holes in global AdS space are reproduced by the equations of non-linear 
fluid dynamics [1] . We refer the reader to [5] for a list of developments and 
references. 



1.1.2 Preview of the fluid/gravity correspondence 

Having provided the reader with a broad, albeit abstract, rationale to as- 
sociate the dynamics of Einstein's equations to that of a quantum field 
theoretic stress tensor, we now provide some specifics that set the stage for 
our discussion. 

According to the gauge/gravity dictionary, distinct asymptotically AdS 
bulk geometries correspond to distinct states in the boundary gauge theory. 
The pure AdS geometry, i.e., the maximally symmetric negatively curved 
spacetime, corresponds to the vacuum state of the gauge theory. A largqj 
Schwarzschild-AdS black hole corresponds to a thermal density matrix in 
the gauge theory. This can be easily conceptualized in terms of the late- 
time configuration a generic state evolves to: in the bulk, the combined ef- 
fect of gravity and negative curvature tends to make a generic configuration 
collapse to form a black hole which settles down to the Schwarzschild-AdS 
geometry, while in the field theory, a generic excitation will eventually ther- 
malize. Note that although the underlying theory is supersymmetric, the 
correspondence applies robustly to non-supersymmetric states such as the 
black holes mentioned above. In this sense, supersymmetry is not needed for 
the correspondence. 

^ Recall that AdS is a space of constant negative curvature, which introduces a length scale, 
called the AdS scale -RAdSi corresponding to the radius of curvature. The black hole size is 
then measured in terms of this AdS scale; large black holes have horizon radius r+ > HAdS- 
We will be focus on the large black hole limit r^ 2> ^AdSi ^-iid therefore consider the planar 
Schwarzschild-AdS black holes. 



1.1 Introduction 5 

On the boundary, the essential physical properties of the gauge theory 
state (such as local energy density, pressure, temperature, entropy current, 
etc.) are captured by the expectation value of the boundary stress tensor, 
which in the bulk is related to normalizable metric perturbations about a 
given state. It can be extracted via a well-defined Brown- York type proce- 



dure [6] as we review later (see (1.39)). 



At the risk of being repetitive we urge the reader to note the distinc- 
tion between the two separate stress tensors that will enter our analysis. 
In our framework, the hulk stress tensor appearing on the r.h.s. of the bulk 
Einstein's equation is zero if we are only interested in the universal sub- 
sector discussed above. On the other hand, the boundary stress tensor T 
is non-zero; it is a measure of the normalizable fall off of the bulk metric 
at the boundary. Note that the boundary stress tensor does not curve the 
boundary spacetime a la Einstein's equations since the boundary metric ^ab 
is non-dynamical and fixed. We will discuss generalizations that allow for 



non-trivial bulk matter in ^1.6 when we move outside the universal stress 
tensor sector. 

To describe gravity duals of fluid flows, a useful starting point is the map 
between the boundary and bulk dynamics in global thermal equilibrium. In 
the field theory, one characterizes thermal equilibrium by a choice of static 
frame and a temperature field. On the gravity side, the natural candidates 
to characterize the equilibrium solution are static (or more generally station- 
ary) black hole spacetimes, as can be seen by demanding regular solutions 
with periodic Euclidean time circle. The temperature of the fluid is given by 
the Hawking temperature of the black hole, while the fluid dynamical veloc- 
ity is captured by the horizon boost velocity of the black hole. For planar 
Schwarzschild-AdS black holes the temperature grows linearly with horizon 
size; the AdS asymptotics thus ensures thermodynamic stability as well as 
providing a natural long wavelength regime. 

Now let us try to gently move away from the equilibrium configura- 
tion. Starting with the stationary black hole (namely the boosted planar 
Schwarzschild-AdSrf+i) solution, we wish to use it to build solutions where 
the fluid dynamical temperature and velocity are slowly-varying functions 
of the boundary directions. Intuitively, this mimics patching together pieces 
of black holes with slightly different temperatures and boosts in a smooth 



way so as to get a regular solution of (1.1). In order to obtain a true so- 
lution of Einstein's equations, the patching up procedure cannot be done 
arbitrarily; one is required at the leading order to constrain the velocity and 
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temperature fields to obey the equations of ideal fluid dynamicsjj Further, 
the solution itself is corrected order by order in a derivative expansion, a 
process that likewise corrects the fluid equations. All these steps may be 
implementeq^ in detail in a systematic boundary gradient expansion. The 
final output is a map between solutions to negative cosmological constant 
gravity and the equations of fluid dynamics in one lower dimension, i.e. the 
fluid/gravity map. 

A noteworthy aspect of this construction is that Einstein's equations 
become tractable due to the long wavelength regime without losing non- 
linearity. From the boundary standpoint, one encounters domains of nearly 
constant fluid variables; these domains can then be extended radially from 
the boundary into the bulk and in each such bulk 'tube', illustrated in Figure 
|1.1[ we are guaranteed to have a solution which is close to the equilibrium 
form. Lest the reader be led astray, we should note that the solutions we 
construct are perturbative and hence approximate. Nevertheless, they are 
'generic' slowly-varying asymptotically-AdS black hole geometries, with no 
Killing fields. 





Figure 1.1 Penrose diagram of the uniform planar black hole (1.17) and the 



causal structure of the spacetimes dual to fluid dynamics illustrating the 
tube structure. Dashed line in the second figure denotes the future event 
horizon, while the shaded tube indicates the region of spacetime over which 
the solution is well approximated by a corresponding uniform black hole. 



Thes e constraints are actually the radial momentum constraints for gravity in AdS and imply 
| |1.2| l. In contrast to the conventional ADM decomposition, we imagine foliating the spacetime 
with timelike leaves and 'evolve' into the AdS bulk radially. 

In the technical implementation of this program, it is important that one respect boundary 
conditions. We require that the bulk metric asymptote to -jab (up to a conformal factor) and 
further be manifestly regular in the part of the spacetime outside of any event horizon. 
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A remarkable outcome of the association between generic black holes and 
fluid flows is that it automatically provides a sensible entropy current with 
non-negative divergence for hydrodynamics. On the gravitational side, en- 
tropy is naturally associated to the area of the event horizon; by pulling 
back this area form to the boundary, we can equip our fluid with a canoni- 
cal entropy current. 

We now turn to the technical aspects of the fluid/gravity map. Following 
a review of fluid dynamics and the perturbative construction of gravitational 
solutions, we finally present the main results (in particular the bulk metric 
and the boundary stress tensor, to second order in boundary derivatives) in 



{lA. The subsequent sections are devoted to describing some implications 



and extensions of the basic construction. 



1.2 Relativistic fluid dynamics 

To set the stage, let us start by reviewing fluid dynamics, explicating the use 
of gradient expansion as an organizational principle. At high temperatures 
every non-trivial quantum field theory (and every experimentally realizable 
system) equilibrates into a fluid phase, i.e., a translationally invariant phase 
in which adiabatic displacement of neighboring elements requires no force. 
Weakly interacting fluids are composed of a collection of a large num- 
ber of long lived partonic excitations which continually collide with each 
other. The time and space intervals between successive collisions of a given 
parton are called the mean free time, t^, and the mean free length, i^, 
respectivelyjj Such fluids are characterized by a parton density function in 
phase space, and the time evolution of this function is governed by the well 
known Boltzmann transport equations of statistical physics. These equa- 
tions have an interesting property: arbitrary initial density functions relax 
to local thermal equilibrium over a time scale of order the mean free time. In 
other words, for t ^ tm, the parton distribution in momentum space approx- 
imately reduces, at every point x, to an equilibrium distribution. However 
the parameterqj characterizing this equilibrium configuration, the tempera- 
ture field T{x) and fluid velocity field u"(x), vary on a length scale that is 
large compared to i^n- T{x) and u"'{x) are the effective dynamical variables 
of the system at later times; their evolution as a function of time is governed 
by the equations of fluid dynamics. 

^ In a relativistic system of massless particles, like Af = 4 SYM, t^ ~ dm- We will assume this 

is the case in our discussions below. 
^ For simplicity, in this discussion we assume that the system has no conserved charges other 

than the stress-energy-momentum tensor, and no other Goldstonc-like light degrees of 

freedom. We discuss generalizations below. 
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Now it turns out that the equations of fluid dynamics may also be derived 
in a much simpler and more general manner, and so apply even at strong 
coupling. The main assumption that underlies fluid dynamics is that systems 
always equilibrate locally over a finite time scale that we continue to refer to 
as im- While this assumption is true of the Boltzmann transport equations, it 
is believed to hold more generally also for strongly coupled fluids. It follows 
immediately that T{x) and u"'{x) are the effective variables for dynamics at 
length and time scales large compared to i^ and tm- As we will now see, the 
equations of fluid dynamics follow inevitably out of this conclusion. 



1.2.1 The equations of fluid dynamics and constitutive relations 

The stress tensor in any d-dimensional quantum field theory on a background 
with metric ^ab obeys the d conservation equations 

VaT^^ = 0. (1.3) 

These equations do not constitute a well defined initial value problem for the 
stress tensor in general as, in d > 2, we have more variables (the ^ d{d + 1) 
independent components of the stress tensor) than equations Fl In the fluid 
dynamical limit, however, the stress tensor is determined as a function of d 



variables, T{x) and n"(x). Consequently, (1.3), supplemented with a formula 
for T as a function of thermodynamical fields, constitute a complete set 
of dynamical equations. These are the equations of fiuid dynamics. 

A constitutive relation that expresses T as a function of T(x), n"(x), and 



their derivatives turns (1.3) into a concrete set of fluid dynamical equations. 
In thermal equilibrium the stress tensor T""^ is given by 

T''^ = {P + p) u" u^ + P -f''^ (1.4) 

where P is the pressure of the fluid and p is its energy density. Recall that 
both p and P are known functions of temperature (determined by the ther- 
modynamic equation of state of the fluid) . For a fluid in local thermal equi- 



librium, (1.4) generalizes to 



T''\x) = [P{x) + p{x)] u"(x) u\x) + P{x) 7"'' + n"^(x) (1.5) 

where P{x) = P{T{x)), p{x) = p{T{x)) and 11"^ (x) represents the contribu- 
tions of derivatives of T{x) and u"'{x) to the stress tensor. This dissipative 

^ For the special case of confomial field theories the number of variables is reduced by one, as 
the trace of the stress tensor vanishes. In this case, at d = 2, we have as many variables as 
equations. This observation underlies the special simplicity of CFTs in d = 2. 
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part n" may then be expanded as 



n'^^ = X;cng) (1.6) 

n=l 

where 11?^-, is defined to be of n order in derivatives of the fluid dynamical 
fields. Note that magnitude of n?''^ relative to the ideal fluid stress tensor is 

b (n) 

approximately (^m/-^)" where L characterizes the length scale of variation 
of the temperature and velocity fields; consequently terms at higher values 
of n are increasingly subdominant in the fluid dynamical limit. 

The explicit form of the functions n?\ can only be derived from a detailed 
study of the dynamics of the specific system. However the allowed forms for 
constitutive relations are significantly constrained by symmetry and other 
general considerations. At first order, for instance, it is possible to assert on 
very general grounds that 



TT(afc> _ pa pb TTcd pab p rrcd _ 9 ab 

(1) = c-r dii(i) - ^^— — j-^ r-cdi-^i) - -^-qa 
1 BP 

j—^Ii''^)Pab-^^{UaU,Tl'^^^) = -Ce 



(1.7) 



where 



p<^b ^ ^a^b ^ ^ab (^_g) 

is the projector onto space in the local fluid rest frame, 

^ab ^ ^{a ^b) _ pac pbd Uj^^ ^^^ _ _1_ p^^ g\ (^^g) 

is the fluid shear tensor, 6 = 'Vcu'^ is the expansion, and the brackets around 
the indices {ab) denote the symmetric transverse traceless part of the ex- 
pression. Here r] and C are arbitrary functions of the temperature, referred 
to as the shear and bulk viscosity, respectively. 



The equation (1.7) is a physically complete specification of the constitu- 



tive relations at first order, even though it leaves P" ^^^ u^ and one linear 
combination of H^-^ P^i, and Ua u^ n?^^ unspecified. This is because T{x) and 
u'^{x) have no intrinsic definition out of equilibrium. All equations of fluid 
dynamics must be 'field redefinition invariant' (invariant under redefinitions 
of T{x) and n"(x) that reduce to identity in equilibrium), and it turns out 



that the l.h.s. of (1.7) are the only field redefinition invariant data in H^n. 



(!)■ 
ab 



The other components of H^^n can be modified at will by an appropriate 
field redefinition, and have no physical significance. 
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The r.h.s. of the two equations in (1.7) represent the most general in- 
equivalent 'tensor' and 'scalar' data that can be constructed from a single 
derivative of fluid dynamical fields compatible with the conservation equa- 
tion at first order which imply that u^'VaT oc 9. 

It is sometimes convenient to fix the field redefinition ambiguity by giving 
the fields u°'{x) and T{x) unambiguous (but arbitrary) meaning. In the so- 
called 'Landau Frame' this is achieved by asserting that, at each point, 

Tl{x) u\x) = -p{x) u'^ix) . (1.10) 

This relation defines ti° by identifying it with the unique timelike eigenvector 
of the stress tensor at any point, and defines the temperature by identifying 
the corresponding eigenvalue with the energy densityj^In the Landau frame, 
which we adopt for most of this chapter, ( |1.7[ ) simplify to 

Il''^^ = -2 7]a''^ -COP"^ . (1.11) 

As the equations of fiuid dynamics are both local and thermodynamical 
in nature, they must respect a local form of the second law of thermody- 
namics. It follows that the equations of fiuid dynamics must be accompanied 
by an entropy current whose divergence is pointwise non-negative in every 
conceivable fluid flow. At first order, for a charge-free fluid, the constraints 
imposed by this requirement are a relatively mild set of inequalities on ry 
and Q: It turns out that the entropy current is constrained to take the form 

J- = su--^u,Iill), (1.12) 

where s is the entropy density. It is possible to demonstrate that the current 



in (1.12) is field redefinition invariant to first order. Note that the second 



term on the r.h.s. of (1.12) vanishes in the Landau frame. Using the Euler 
relation (p+P = sT) and the Gibbs-Duhem relation (dP = sdT) of thermo- 
dynamics along with the equations of motion, it follows that the divergence 
of this entropy current is given by 

VaJ! = -Va (I) n^i') . (1.13) 



Using (1.7) (or more simply (1.11) in the Landau frame), it is easy to verify 
that positivity of the entropy current requires that r/ > and C ^ 0. At 
higher orders in the derivative expansion (for uncharged fluids) and even at 
first order for more complicated fluids (e.g. charged fluids and superfluids) 
the requirement of positivity of the entropy current imposes more than a 

* Note that the equations l|1.10[| are true in equilibrium. 
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set of inequalities on transport coefficients; it forces Unear combinations of 
otherwise arbitrary transport coefficients to vanish [7]. 

In the rest of this chapter we will be especially interested in the fluid 
dynamics of conformal field theories. These theories enjoy three key simpli- 
fications. First, as they have no dimensionless parameters, the dependence of 
all physical quantities (e.g. P, p, rj) on temperature follows on dimensional 
grounds. In particular, 

P = aT'^, p={d-l)aT'^, rj = rj' T'^-^ (1.14) 

where a and rj' are dimensionless constants. Second, the stress tensor in any 
CFT is necessarily traceless. It in particular follows from this condition that 
C = 0. Finally, the stress tensor in such theories must transform covariantly 
under Weyl transformations. This imposes additional restrictions on the 
stress tensor at higher orders in the derivative expansion. In summary, for 
a conformal fluid, the stress tensor up to first order has the form 

j.ab ^^rj^d fdu^u'' + -f'''') - ri'T'^-^a'''' (1.15) 



where a and r/' are pure numbers and ry' > 0. 

The constraints on allowed forms of the constitutive relations at higher 
order are more complicated. The most general allowed equations of second 
order fiuid dynamics have largely (but perhaps not completely) been worked 
out in [8] . In the following we will determine the second order fiuid equations 
for A/" = 4 Yang Mills at strong coupling using the fiuid/gravity duality. 



1.2.2 The Navier- Stokes scaling limit 

An interesting fact about the equations of relativistic (or any other compress- 
ible) fiuid dynamics is that they reduce to a universal form under a combined 
low amplitude and long wavelength scaling. Consider a uniform fiuid at rest, 
perturbed so that the amplitude in velocity fiuctuations is small (scales like 
e) and the amplitude in temperature fiuctuations is smaller (scales like e^). 
We also require that the wavelength of spatial fluctuations is large (scales 
like 1/e) and that their temporal scale even larger (scales like 1/e^). We then 
take the strict e — t- limit. In this limit: 

1) The fluid is non- relativistic, as all velocities are parametrically smaller 
than the speed of light. 

2) The fluid is incompressible, as all velocities are parametrically smaller 
than the speed of sound (recall that a sound wave is a compression wave, 
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and that fluid flows at velocities smaller than the speed of sound are 
effectively incompressible) . 

3) The temporal component of the energy conservation equations reduces, 
at leading order ©(e^) to the continuity equation V ■ v = 0. We use the 
symbol v for the non-relativistic spatial velocity. 

4) The spatial component of the energy conservation equations reduce at 
leading order, 0(6^), to the famous (non-relativistic) Navier-Stokes equa- 
tions 

ij+v-Vv = -VP + uV'^v (1.16) 

with kinematic viscosity 

Po + Po' 

where po and Pq are the background values of the density and pressure of 
the fluid. 

Note that the Navier-Stokes equations are homogeneous neither in ampli- 
tude of fluctuations (the convective term is non-linear), nor in derivatives 



(the viscous term is quadratic in derivatives). All terms retained in (1.16) 
are equally important in the e — )• limit; in particular, the parameter ly may 
be set to unity by a uniform rescaling of space and time. In contrast, the 
viscosities ry, C give a subleading correction to ideal fluid dynamics in the 
derivative expansion of relativistic fluid dynamics. 



Taking the spatial divergence of (1.16) one sees that the pressure may 



be solved for in terms of the velocity field on any given time slice; so the 
pressure is not an independent degree of freedom. The initial data of the 
Navier-Stokes equations comprise just the components of a divergence-free 
velocity field specified on any time slice. This reduction is not surprising 
since the sound waves are being projected out in this limit. 



The incompressible Navier-Stokes equations (1.16) describe a wide variety 
of sometimes extremely complicated phenomena like turbulence. Despite the 
fact that these equations have been studied for almost 200 years, their non- 
linear phenomenology remains rather poorly understood. One of the hopes 
of the fluid/gravity correspondence is to shed a new (geometric) light on 
some of these issues. 



1.3 Perturbative construction of gravity solutions 

We have seen that fluid dynamics can be treated systematically as the theory 
of long wavelength fluctuations about thermal equilibrium. We are now going 
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to construct gravitational solutions dual to fluid flows by formalizing this 
intuition to set up an algorithmic procedure to construct slowly varying 



dynamical black hole spacetimes as solutions to (1.1) 



1.3.1 Global thermal equilibrium from gravity 

The starting point from the gravitational perspective is a solution that corre- 
sponds to global thermal equilibrium. For the moment let us consider a con- 
formal field theory on Minkowski space jab = Vab- From the gauge/gravity 
correspondence we know that the dual geometry in the bulk is the planar 
Schwarzschild-AdSd+i geometry, whose metric in static coordinates (with 
x"" = {t, y*}) is given by 



dr'^ 9 . , ,■ , ,• „, ^ -, /47r 



- 



(1.17) 
This is a one-parameter family of solutions, parameterized in terms of the 
black hole temperature T, which determines the horizon radius, r+ = ^g— , 
where / vanishes. It is easy to generate a d-parameter family of solutions 



by boosting (1.17) along the translationally invariant spatial directions y', 
leading to a solution parameterized by a (normalized) timelike velocity field 
u"". The parameters which characterize the bulk solution are precisely the 
basic hydrodynamical degrees of freedom, viz., temperature T and veloc- 
ity u°' of the black hole. It is easy to see that the solution induces on the 
Minkowski boundary of the AdS^+i spacetime a stress tensor which precisely 



takes the ideal fiuid form of (1.4) with thermodynamic parameters specified 
by (1.14). The normalization constant a is fixed by the gravitational theory 



to be a = TrfTTT- 

Consider this d-parameter family of boosted planar Schwarzschild-AdS^+i 
geometries, each of which holographically map to an ideal (conformal) fiuid 
living on R '^ (endowed with the Minkowski metric). This fiuid is in 
global thermal equilibrium and one should be able to describe the long 
wavelength but arbitrary-amplitude fluctuations away from equilibrium via 
hydrodynamics. This class of fluctuations has a bulk geometric avatar; we 
now describe an algorithmic procedure that enables us to construct such 
asymptotically locally AdS^+i black hole geometries which are generically 
inhomogeneous and dynamical. This procedure crucially relies on the fact 
that hydrodynamics, as indeed any effective field theory, can be systemati- 
cally studied in a gradient expansion. 
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1.3.2 The perturbation theory 

We start by considering perturbations to the seed geometry characteriz- 
ing equihbrium. Take the boosted planar Schwarzschild-AdS^+i spacetime 



(1.17), for convenience rewritten in ingoing coordinates so as to remove the 
coordinate singularity on the horizon, and replace the parameters Ua and T 
by functions of the boundary coordinates x", 

ds^ = -2 Ua{x) dx" dr-r'^ f {r/T{x)) Ua{x) Ub{x) dx"- dx^+r'^ Pab{x) dx" dx^ , 

(1.18) 



with /(r) as specified in (1.17) and Pab given by (1.8) with ^ab = ilab- This 
metric, which we henceforth denote as g\j,u{T{x),u°-{x)), is not a solution to 
Einstein's equations. It however has two felicitous features: (i) it is regular 
for all r > and (ii) if the functions T{x) and u'^{x) are chosen so as to 
have small derivatives, then it can be approximated in local domains by a 
corresponding boosted black hole solution. These observations lead us to 



consider an iterative procedure to correct (1.18) order by order in a gradient 
expansion. We will find that we however cannot specify just any slowly 
varying T{x) and u"'{x) (recall that x"" included the temporal direction). A 
true solution to Einstein's equations is obtained only when the functions 
T{x) and u°'{x) in addition to being slowly varying satisfy a set of equations 
which happen to be precisely the conservation equations of fluid dynamics. 
Let us record that we have fixed gauge by setting g^r = and Qra = —Ua- 

Since we want to keep track of the derivatives with respect to the boundary 
coordinates, it is useful to introduce a book-keeping parameter e and regard 
the variables of the problem as functions of rescaled boundary coordinates 
ex"". At the end of the day e may be set to unity. With this in mind, let us 
consider the corrections to the seed metric in a gradient expansion: 



9^.^ = E e''9^^KT{ex),uyEx)) , u'^ = E^'^^^'^M ^ ^ = E^'^^'^(^^) 

fc=0 fc=0 fc=0 

(1.19) 
where the correction pieces g\_iJ , u""^^' and T^^' are to be determined by 
solving Einstein's equations to the A;*^ order in the gradient expansion. The 



ansatz (1.19) should therefore be inserted into the Einstein's equations (1.1) 
and the result expanded in powers of e. 

Let us examine the resulting structure in abstraction first. For the sake of 
argument, assume that we have determined gj^ for m <n — l and T^^' and 
ya(m) £qj, 777, < ri — 2. At Order e" one finds that Einstein's equations reduce 
to a set of inhomogeneous linear differential equations whose structure can 
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be schematically written as 
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(o)(rW,n 



a(0)^ 



9 



(n) 



(1.20) 



where we have dropped the spacetime indices for notational clarity (c.f., 
[U [5] for the explicit equations). Since each derivative with respect to x" 
is accompanied by a power of e, it follows that the linear operator H is 
constructed purely from the data of the equilibrium Schwarzschild-AdSd+i 
geometry. This means that HI is at most a second-order differential operator 
with respect to the radial variable r. Moreover, it has to be independent of n. 
Thus the perturbation theory in e is ultra-local in the boundary coordinates, 
implying that we can solve the equations of motion of the bulk spacetime 
point by point on the boundary! 

On the right hand side of (1.20) we collect all order e" terms which do 
not have explicit radial derivatives into a source term s^, which is then a 
complicated construct involving contributions from different orders in per- 
turbation theory. It is a local expression of (n — m) order in boundary 
derivatives of T^'^' and u"'^^) for m <n — l, and ascertaining it is the most 
substantial part of the computation. 



The reader may be puzzled by the following aspect of (1.20): while we 



d{d+l) 



d(d-l) 



have -^ — - equations, we have only ^ ^ ' variables after fixing the gauge 
redundancy. This implies that a subset of Einstein's equations has a distin- 
guished status as constraint equations, while the remainder are the physical 
dynamical equations. 



To understand this let us examine the differential equations (1.20) by in- 
voking the canonical split of our bulk coordinates X^ = (r, x"). The Era 
equations are the momentum constraint equations for 'evolution' in the ra- 
dial direction. These equations are special in several ways. To start with, 
they need only be satisfied on a single r slice; the 'dynamical' equations 
{Eab) then ensure that they will be solved on every r slice. For this reason, 
it is consistent to study these equations just at the boundary, where they 
turn out to reduce merely to the equations of conservation of the boundary 
stress tensor 

VaT^ti) = • (1-21) 

(See Q.4.4 for the definition of the boundary stress tensor.) Note that at 



n order the equations (1.21) depend only on the boundary stress tensor 
built out of the spacetime metric at order n — 1. This is because (1.21) 



has an explicit boundary derivative which carries its own effective power 
of £. The net upshot is that the unknown metric g^"^' does not enter the 



equations ( 1.21 ) at all (the operator EI in (1.20 ) vanishes for these solutions) 
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Hence, (1.21) is instead a constraint on the solution already obtained at 
one lower order in perturbation theory. As we will see below, the solution 
for ^'"■) of the dynamical equations at each order in perturbation theory is 
uniquely obtained in terms of the previous order solution, and so, ultimately, 
in terms of the velocity and temperature fields that enter the starting ansatz 



(the zeroth-order term in perturbation theory). Consequently, (1.21) is an 
equation which constrains the starting velocity and temperature fields, and 
turns out to be the equation of boundary fluid dynamics. 

The remaining equations E^r (the 'Hamiltonian constraint' for radial evo- 
lution) and Eab are dynamical equations with the operator H being a second 
order differential operator in r. Exploiting the spatial rotational symmetry 
of the seed solution, these equations can be decoupled and solved by quadra- 
tures, 

g^"> = particular(s„) + homogeneous(]HI) . (1-22) 

A unique solution to the dynamical equations is obtained upon specification 
of boundary conditions: normalizability at infinity and regularity in the in- 
terior for all r > 0. These turn out to specify the solution completel}|j and 
one ends up with a regular black hole geometry at each given order in the e 
expansion. 

In summary, at any order in the perturbative expansion one solves the 
constraint equations, enforcing fluid dynamical equations on the 'initial' 
data. One then solves for the corrected metric. This correction feeds into 
the constraint equations giving corrected equations of fluid dynamics, and 
so on. The process may be iterated to any desired order, thereby yielding a 
systematic derivative expansion of the equations of fluid dynamics. 



1.4 Results at 2nd order 

Having seen abstractly the iterative procedure which perturbatively corrects 
the seed metric to obtain a solution to Einstein's equations at arbitrary order 
in the gradient expansion, we now turn to the results of this construction (for 
now still considering only energy- momentum transport on the boundary). 
While our discussion so far has been restricted to the case of a flat boundary 
metric 'jab = ijab^ the observation we made about the ultra- locality of the 
perturbation theory allows us to immediately generalize to slowly-varying 
curved boundary metrics. Given a metric 'jab on the boundary, we can ex- 
ploit the freedom to pass over to a Gaussian normal coordinate chart about 

^ Modulo the fact that the operator H has zero modes which are to be accounted for by 
re-definitions of the background values of T and m". 
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the point under consideration, and account for the curvatures which arise 
starting with the second order in the e expansion via the computation of 
appropriate source terms. We will therefore present the results below for this 
more general setting. Before we do so, we will take the opportunity to review 
a beautiful technical framework developed by [S] to simplify the results for 
conformal fluids. 



1.4-1 Weyl-covariant formalism 

The vacuum AdS^+i spacetime is dual to the vacuum state of a conformal 
field theory. If we are interested in the hydrodynamic description of the latter 
on a background manifold Bd, then rather than focusing on the metric 7^^ 
of this geometry, we can consider the conformal class of metrics {Bd, C). On 
this conformal class there is a natural derivative operator, defined through 
a Weyl connection, which efficiently keeps track of Weyl transformation 
properties of various operators. This is all the more natural in the context 
of fluid dynamics where there is a distinguished vector field, the velocity u"", 
defined to be the (normalized) timelike eigenvector of the stress tensor. 

Let us first start with local Weyl rescalings of the boundary metric which 
transforms homogeneously, i.e.. 

We will call a tensor Q with components Qa^.'.'.'a^ conformally covariant and 
of weight w if it transforms ho moge neously under Weyl rescalings of the 



metric, i.e., Q = e~'^'^ Q under (1.23). The velocity field u"" transforms as a 
weight 1 tensor while the stress tensor T of a conformal fluid has weight 
{d + 2) in d-spacetime dimensions. 

One defines a class of torsionless connections, called the Weyl connec- 
tions, characterized by a connection one form Aa^ whose associated covariant 
derivative V ^^ captures the fact that the metric transforms homogeneously 
under conformal transformations (with weight —2). In particular, for every 
metric in the conformal class C, 

VryV = 2A76c. (1.24) 

Given this derivative structure, we can go ahead and define a Weyl co- 
variant derivative Da = Va '^^ + w Aa which is metric compatible and 
whose action on tensors transforming homogeneously with weight w (i.e., 
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Qa-- =£-«-</. ga.-) is given by 



jcb A" - 5i A, - 5t Ac) Qf. -■■■ . (1.25) 



The connection has been defined so that the Weyl covariant derivative of 
a conformally covariant tensor transforms homogeneously with the same 
weight as the tensor itself. 

In hydrodynamics we will require that the Weyl covariant derivative of 
the fluid velocity be transverse and traceless, 

t("Paii^ = 0, Vau"' = Q, (1.26) 

which enables one to uniquely determine the connection one-form Aa to be 
the distinguished vector field 

Aa = u"" Vc Ua - J— Ua Vc u" = Qa - j— Ua , (1.27) 

built from the fluid velocity field. 

One can rewrite the various quantities appearing in the gradient expansion 
of the stress tensor in this Weyl covariant notation. For instance, at first 
order in derivatives, we have the shear and vorticity constructed from the 
velocity field: 

^ab ^ jy(a^b) ^ ^ab ^ _p[a^6] ^ (^ 38) 

both of which have weight w = 3. The fluid dynamical equations, viz., 
stress tensor conservation, are simply P^ T""^ = in this Weyl covariant 



language (which is equivalent to (1.3) since (1.25) with w = d + 2 gives 



Va T"-^ = Va T"-'^ + T^"- A^ and the conformal fluid stress tensor must be 
traceless) . 



1.4-2 Generic asymptotically AdS black hole metric 

We now have at our disposal all the technical machinery necessary to present 
the results for the gravity dual of non- linear fluid dynamics. By a suitable 
choice of gauge (a slight generalization of the Eddington-Finkelstein coordi- 
nates), one can express the bulk metric g^y in the form 

ds^ = -2 Ua{x) dx" (dr + 5Jfe(r, x) dx^) + (5ab{r, x) dx" dx^ , (1.29) 
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where the fields QJa and <3ab are functions of r and x° which admit an 
expansion in the boundary derivatives. In the parameterization used in [10] 
one finds the metric functions are given up to second order in derivatives as: 



^, = rAa- Sac u' - Mr/T) P: V^ a\ 



\ r^ f{r/T) + i (1 - fir/T)) c.^.u^'^ + ^2{r/T) ^ 



cd 



<5ab = r" Pah - < ^cb + 2 {r/Tf gi (6 r) 

cd 

-92{br)^^Pab-Q3{r/T) 

+ 94{r/T) [llab + '^iab]- 



AttT 



crab + Qiir/T)a^''a, 



cb 



^lab + 7^ 2^3a6 + 2 ^2ab 



(1.30) 

Here 5a;, = j^^g ( ^afe — 2{d-i) ^"^ ) ^^ ^^^ Schouten tensor of the boundary 
metric, where the Weyl covariant curvature tensors are 



T^ah = Rab -{d-2) {Va A + A A " A^ lab) " Qab^c A - Tab ,^ ^^. 

7^ = 7^a" = i? - 2 (d - 1) Vc A + (d - 2)(d - 1) A 



with Fab = Va ^5 — Vb Aa- Apart from the shear and vorticity tensors ( 1.28 ) 



constructed from the fluid velocity, we also encounter four of the five second 
order tensors which form a Weyl covariant basis. 






crab 

crab _ 
-^4 - 



C 



acbd 



UcUd , 



2a''^''uj^l, %^ = >A<',.,b) 



:i.32) 



IjJ ^ UJ 



Note that the tensor &ab is clearly transverse, since it is built out of op- 
erators that are orthogonal to the velocity, and it can be inverted via the 
relation (©~^) ®c6 = P^ ■ The induced metric on the boundary in these 
coordinates takes the form: 



7a6 = lim -^ {(5ab-2u,a^i,)) 



(1.33) 



which is crucially used to raise and lower the boundary indices. 

Finally, the various functions gj and Dj appearing in the metric are given 
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in terms of definite integrals once one has inverted tlie operator H: 



01 (y) 



c 



d-1 



1 



, '^'^C(C'^-i 



92iy) = 2y' f^f dCC's'iiO 



5 

^d-2 _ 1 






fi(y) 
02(y) 






1 l°° di 



i-e(e-i)0;(e)-2(d-i)^'^-i 



+ (2(d-i)e"-(d-2)) I'dcc'sUc) 



(1.34) 
The asymptotic behavior of these functions Qi{r/T) and Vi{r/T) is important 



for the stress tensor computation of ^ 1.4.4 and can be found in [lOJ. 



1.4-3 Event horizon and entropy current 



The metric (1.29), (1.30) solves Einstein's equations to second order in the 



gradient expansion, provided the first order stress tensor (which takes the 



form (1.15) with the coefficients extracted from the 1st order bulk metric. 



and given explicitly below in { 1.4.4) satisfies the hydrodynamic conservation 



equations. While this already establishes a firm connection between solutions 
of Einstein's equations and those of fluid dynamics (in one lower dimension), 
it is imperative to establish that the bulk geometry we describe is regular 
everywhere outside the curvature singularity at r = 0. 

Although one can utilize the behavior of the metric functions and iter- 
atively argue that the sources are regular order by order in perturbation 
theory, it is convenient to establish once and for all that what one has con- 
structed is a black hole spacetime with a regular event horizon. Doing so 
involves ascertaining the location of the event horizon. A-priori, this sounds 
like a tall order, especially given that explicit solution is contingent on hav- 
ing solved the fluid equations. Moreover, as is well known, the event horizon 
is a teleological concept (it is the boundary of the past of future null in- 
finity) whose determination requires knowing the entire future history of 
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the spacetime. However, with one key assumption of late-time relaxation 
which is natural from fluid dynamics, it turns out to be possible to deter- 
mine the location of the event horizon locally within our gradient expansion. 
Apart from showing regularity, this has the additional virtue of enabling us 
determine a natural entropy current for fluid dynamics [TT] . 

Since generic flows of dissipative fluids tend to approach global equilib- 
rium at late times, it follows that the corresponding event horizon has to 
approach the radial position determined by the local late-time temperature 
of the fluid. In particular, we look for a null co-dimension one surface given 
by the equation 5-^(r, x) = r — r-^(x) = with the correct asymptotics. 
The function r-^(x) should be parameterized within the gradient expansion 
ry_{x) = ^ J^^' + Y^j. e^r(fc)(x). The corrections r(fc)(x) are determined by 
solving the null condition g^^'^ d^Sy, d^.Sy, = 0. The resulting equations are 
algebraic for r^^) and to second order in gradients one finds that, for the 



solution ( 1.29 )-( 1.34) 



4ttT{x) d 

2(^2 + d- 4) 2 02(1) 



rw(x) = ^^^^ + ^^^^ [^iaatc7-'' + ^2u;abU^''' + ^3n} (1.35) 
with 



«i 



d"^ {d - 1) {d - 2) d{d-l) 



K ^ + 2 1 

^ 2d{d-2) ' ^ d{d-l){d-2) 



(1.36) 



This indeed establishes that the solutions we have constructed in ilA.2 to 
2nd order qualify to be called inhomogeneous, dynamical black holes. 

Note that in general, beyond the leading order, the horizon position and 
generators are not simply given by the corresponding fluid temperature and 
velocity (for example, while the horizon generators must be vorticity-free, 
coab need not vanish for the boundary fluid). In some sense, while the black 
hole horizon is distinguished in the bulk, physics appears simpler when ex- 
pressed in terms of the fluid data living on the boundary. 

Having determined the event horizon of the gravity solution, we immedi- 
ately have access to an important hydrodynamic quantity, viz., the entropy 
current. For a black hole spacetime it is natural to view the area of the event 
horizon as an entropy a la Bekenstein-Hawking |121ll3j. In fact, by suitably 
foliating the event horizon with spatial slices (propagated forward by the 
null generator), we can equivalently talk about an area {d — l)-form a-^ on 
these slices. Since we imagine the dual fluid living on the boundary of the 
spacetime, it is natural to pull-back this area form out to the boundary. A 
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canonical choice is to pull-back along radially ingoing null geodesies jllj . 



which is quite easy to implement for the metric (1.29), where the lines of 
x"- = constant are precisely such geodesies. We then have a {d — l)-form on 
the boundary which can be dualized to a one-form or equivalently a current 
Jg, which is the entropy current on the boundary. Not only does this defini- 
tion agree with the equilibrium notion of entropy of the fluid, but also thanks 
to the area theorem of black hole horizons, we are immediately guaranteed 
that this current has manifest non-negative divergence as demanded by the 
second law. The hydrodynamic entropy current takes the general form 

,2 ^ (1-37) 



+ TTTT^ (^1 ^- ^" + B2 Vc i^n + 



sd 

where s is the entropy density and ^1,2,3 1 ^1,2 are a- priori arbitrary nu- 
merical coefficients. While 'DaJg > only demands that Bi + 2A3 = 0, 



the gravity solution (1.29) fixes all the coefficients in (1.37) explicitly. In 
particular, we obtain 

1 /47rTY"^ , 2 , , , /I , , 2 , / 

Ao = ; , -Bi = — 2 A'i = — r- , Bo = — . 

2d' ^ d{d-2) ' d-2 

(1.38) 

We should note here that there is an ambiguity in pulling back the area- 
form from the event horizon to the boundary, for one can supplement the 
pull-back map with a boundary diffeomorphism, which affects the coefficient 
Ai above. Since this just relabels boundary points in the gradient expansion, 
one is tempted to think of this ambiguity as unphysical. However, it is rather 
curious that if one tries to pull-back the area form from quasi-local horizons 
one encounters a shifted value of Ai [14j , which suggests that there is perhaps 
more to this ambiguity than meets the eye. 



1.4-4 Stress tensor of dissipative fluid 

Given an asymptotically locally AdS^+i metric, one can construct a quasi- 
local boundary tensor which is manifestly conserved and is associated with 
the stress-energy-momentum tensor of the conformal field theory |15| [6] . To 
perform the computation one regulates the bulk spacetime by introducing an 
explicit cut-off at r = Too . The boundary stress tensor is given in terms of the 
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extrinsic curvature Kab of this surface, defined in terms of its unit outward 
pointing normal n" as Kab = lac^^^b- In addition to tlie extrinsic curvature 
one also has contributions from the counter-terms necessary to obtain a finite 
boundary stress tensor. Denoting the curvatures of the boundary metric by 
'^R etc., this is given (to 2nd order) as 



tab 
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(d+i) 
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Kab - K-fab + id- 1) -fab 



m. 



1 



ab 



^ ^Rlab 
(1.39) 

For the gravity duals to fluid dynamics constructed in ^1.4.2 one finds 



that the boundary stress tensor takes the form (1.5) with the dissipative 



part, at the first and second order, given by this gravitational construction 
to be 



jrab 



n' 



ab 

(2) 






(1.40) 



where the tensors Tf^ were defined earlier in (1.32). With the tensor struc- 



ture determined, one is just left with fixing the six transport coefficients, ry, 
Ttt, k, and Aj for i = {1,2,3}, which completely characterizes the flow of a 
non-linear viscous fluid with a gravitational dual. The transport coefficients 
for conformal fluids in d-dimensional boundary turn out to be 
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2-K 



Harmonic 
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T 



2-K {d-2)J 
A3 = 



(1.41) 



where Harmonic (x) is the harmonic number function. Setting d = 4 in the 
above expressions and using the fact that Harmonic(— 2) = —2 log(2) to- 
gether with the replacement — ^^ = ^^ N"^, one can obtain the transport 



coefficients for SU{N) J\f = 4 Super Yang-Mills theory |16tlT]. This has been 
used for real data analysis from e.g. RHIC. 



One immediate consequence of ( 1.41 ) and (1.38 ) is that our fluid saturates 
the famous bound on the viscosity to entropy density ratio, - > j-, |17j . 



This bound is saturated by a large class of two-derivative theories of gravity, 
and it is indeed experimentally satisfied by all presently-known systems in 
nature. Intriguingly, cold atoms at unitarity and quark-gluon plasma both 
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come near to saturating this bound [18]. Its status in more general theories 
is currently under active debate [T9] . 



Moreover, (1.41) reveals further intriguing relations between the coeffi- 
cients, which hint at the specific nature of any conformal fluid which admits 
a gravitational dual. For example, the result that A3 = is universal but 
non-trivial from the fluid standpoint. We also see that 2r]TT^ = 4Ai + A2 
for all d] this in fact was shown to hold quite generally in a large class of 
two-derivative theories of gravity (including matter couplings) |,20j . 



1.5 Specific fluid flovi^s and their gravitational analog 

The construction presented above can be generalized in many interesting 
ways; however before indicating the most important of these in §1.6[ we first 
pause to discuss some of the special cases of the framework explained in 
the preceding section. One of the reasons to discuss these special cases is 
that while we have demonstrated the existence of a map from the equations 
governing fluid dynamics to those governing the dynamics of gravity, we did 
not at any stage solve the fluid equations explicitly. The felicitous feature of 
our construction was the ultra-locality along the boundary directions which 
allowed us to implement the construction in terms of local solutions to the 
conservation equations. Construction of novel fluid flows and generic behav- 
ior of the relativistic conservation equations are interesting (and perhaps 
hard) questions. Nevertheless there are some corners where we can gain an- 
alytic control which serves not only as a check that the fluid/gravity solution 
set is non-empty, but also provides a point of contact with previous studies 
of the hydrodynamic regime in the AdS/CFT literature. 

1.5.1 Linearized setting: quasinormal modes 

Above we have established a map between any solution of the equations of 
fluid dynamics and long wavelength solutions of Einstein gravity with a neg- 
ative cosmological constant. In order to find explicit gravitational solutions 
we need a class of explicit solutions to the equations of fluid dynamics. In 
this subsection and the next we will study such examples. 

It is of course easy to solve the equations of fluid dynamics, derived above, 
when linearized about static equilibrium. Utilizing translational invariance, 
we search for solutions of the form 

with purely spatial velocity fluctuations 6v^ . The resulting linear equations 
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require that the matrix of coefficients M{u},k) annihilate the length-d col- 
umn vector with entries 5v^ and 6T. So the spectrum u}{k) is obtained as 
the roots of the d*^ order polynomial det(M) = 0. At leading (ideal fluid) 
order, the d roots to this equation turn out to be w = ib^T= (the sound 
modes of the fluid) and w = with degeneracy d — 2 (the shear modes of the 
fluid) . These modes and their corresponding eigenvectors receive corrections 
at higher orders in the derivative expansion; in particular the shear modes 
pick up nontrivial k dependence, uj ccik'^, at first order. The explicit solu- 
tions are easily determined (see [IJ). Employing the fluid/gravity map then 



yields explicit linearized solutions of the Einstein equations (1.1) about the 
planar black hole background, whose study in fact predates the fluid/gravity 
map by almost 10 years. 

The spectral problem of linearized fluctuations about a black hole solu- 
tion is of course a well-studied topic (cf., [2T]). It is known that due to the 
presence of a horizon, black holes admit no normal modes. Instead, by im- 
posing regularity at the future event horizon, one finds quasinormal modes, 
i.e., modes which have complex frequencies characterizing decay of pertur- 
bations at late times. Mathematically these are related to the poles of the 
retarded Green's function computed in the black hole background. 

In terms of the gauge/gravity perspective, such quasinormal modes de- 
scribe the timescale for return to thermal equilibrium in the field theory 
|22j . Asymptotically AdS black holes host an infinite family of quasinormal 
modes. All except d of these are 'massive'; their frequency remains finite 
(and has a finite imaginary part or decay rate oc T) even in the limit A; — )• 0. 
However, planar black holes admit exactly d special 'massless' quasinormal 
modes. These so-called hydrodynamic modes can have arbitrarily low fre- 
quency at long spatial wavelengths and therefore fall within the long wave- 
length regime. These are precisely the sound and shear modes described 
above. 

The fact that the dispersion relations of the massless quasinormal modes 
agrees with hydrodynamic dispersion relations was first demonstrated in the 
pioneering works [23l El] , who mapped the Schwarzschild-AdS quasinormal 
modes to sound and shear modes of the dual field theory (perturbed from 
thermal equilibrium). 

1.5.2 Rotating black holes in global AdS space 

A second class of examples which accord analytic control are explicit solu- 
tions corresponding to stationary configurations in hydrodynamics. While 
on fiat space there are no interesting stationary fiows other than a uniformly 
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boosted fluid (whose dual is the seed solution we started with), it turns out 
that by a suitable choice of background geometry one can derive nontriv- 
ial flows. We now describe one such flow on the Einstein Static Universe 
(R X S*^"^) which allows one to make contact with rotating black holes in 
asymptotically (globally) AdS spacetimes. 

Fluid flows of a conformal relativistic d-dimensional fluid on spatial S*^"^ 
conserve angular momentum in addition to energy. The angular momentum 
on S*^"^ is a rank-d antisymmetric matrix, which can be brought to canonical 
form by an SO{d) similarity transform (a rotation) and is therefore labeled 
by its [d/2] inequivalent eigenvalues. For every physically allowed choice of 
these [d/2] angular momenta together with the energy, an arbitrary fluid 
flow eventually settles down into an equilibrium stationary configuration. 

The stationary configurations of viscous conformal fluids on spheres turn 
out to be extremely simple. The velocity field is simply that of a rigid rota- 
tion. Focusing on the case of d = 2n for concreteness, the metric of a unit 
g2n-i jj^^y ]-,g written in terms of the direction cosines /Uj as 

n n 

ds^2n-i = y^fJ-i d(j)i + d^^ , where /,/^i = 1- (l-^^) 

In these coordinates the velocity and temperature fields of stationary flows 
take the form 



(1.43) 
This flow is Weyl equivalent to a uniform velocity and temperature config- 
uration on a confomally rescaled spacetime PP, i.e., 

n 

ds^ = 72 {-dt" + dsl,^.,) , u^^da = dt + Y^ cJid^^ , Ts = To . (1.44) 



It is not difficult to verify that (1.44) provides a non-dissipative solution 



to the equations of fluid dynamics at least to second order (the solution is 
non-dissipative because it is stationary; equivalently the divergence of the 
entropy current vanishes). The constant parameters cjj and Tq of this solution 
turn out to have thermodynamical significance: they are simply the angular 
velocity (chemical potential for angular momentum) and temperature of the 
fluid configuration. 

According to the AdS/CFT correspondence, the dual description of a con- 
formal field theory on MxS'^^-'^ is simply asymptotically global AdS^^+i space. 
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If we pump a large amount of energy and angular momentum into global 
AdS^+i space and let the system relax, we expect the eventual equilibrium 
configuration to be that of a large rotating black hole. This reasoning leads 
to a prediction: the fluid/gravity map applied to ( |1.44[ ) should produce the 
(independently known) metric of large rotating AdS^^+i black hole, expanded 
to second order in the derivative expansion. This predictioq^^ | has been veri- 
fied in detail in the following manner. It is possible to transform the exactly 
known metric of the rotating AdS^+i black holes to the fluid/gravity gauge 



described in Q.3 above. This maneuver in fact turns out to greatly simplify 



the rotating black hole metric which takes the form (|1.29) with 



^a{r,x)=rAa-SacU'' 



2 d? det [r 5\ - uj\\ ' 
*5afe(r, x) = r^ Pab - u}^"" Ucb ■ (1-45) 

The above is an exact rewriting of the rotating AdS black holes (in even di- 



mensions) of [26]. The metric (1.45) may be expanded in derivatives simply 
by expanding the inverse determinant in powers of (jJ. By truncating this 
expansion at second order we recover exactly the metric dual to the fluid 



flow (1.44). It is rather remarkable the full black hole solution can be writ- 



ten economically within the fluid/gravity metric ansatz, perhaps suggesting 



greater utility for metrics of the form (1.29) 



1.5.3 Non-relativistic fluids 

While relativistic fluids are interesting in astrophysical or high energy plasma 
physics contexts, most fluids we encounter in everyday situations are non- 
relativistic. Furthermore, for many practical applications one is usually in- 
terested in their dynamics in the incompressible regime, which is attained 
by projecting out the sound mode. It is natural to ask whether this regime 
is accessible to fluid/gravity; an affirmative answer is suggested by the dis- 
cussion in §1.2.2 namely one only needs to implement the Navier-Stokes 



scaling limit directly in the fluid/gravity solutions. This procedure has been 
carried out in \T7] to obtain the gravitational dual of non-relativistic incom- 
pressible fluid flows. In principle this provides a geometric window to explore 
phenomenologically interesting fluid flows. 

^'^ The first observation that the properties of large rotating blacfc holes should be reproduced 
from fluid dynamics was made in [4], a precursor to the fluid/gravity map. Rotating black 
holes were further studied in 1251 in 4 dimensions, and fully analyzed in d dimensions in [10) . 
Here for illustration we give the general result of ,10_, who find the explicit coordinate 
transformation to rewrite the rotating black hole solution in AdS^+i of 1261 in fluid variables. 
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1.6 Extensions beyond conformal fluids 

The fluid/gravity correspondence was originally derived for the case of con- 
formal fluids which are related to gravitational dynamics in an asymp- 
totically AdS spacetime. Conformal theories are rather special and 'ordi- 
nary' fluids one encounters everyday deviate significantly from this behav- 
ior. Hence it would be useful to look for extensions of the basic framework 
which allows for these generalizations. As already indicated earlier, this can 
be done at the expense of complicating the system of gravitational equations, 
and accompanying loss of universality. Nevertheless, forays into these areas 
have revealed very interesting lessons about fluid dynamics in general that 
transcend the fluid/gravity correspondence. In this section we take stock of 
some of these developments. 



1.6.1 Non- conformal fluids 

The first generalization we consider is provided by a handy trick to obtain 
a particular class of non-conformal theories. It turns out that by exploiting 
the gauge/gravity duality for a special class of theories, viz., theories that 
naturally arise on the world-volume of Dp-branes, one finds a surprisingly 
tractable class of non-conformal fluids. 

Let us first consider the special case of the D4-brane which is a solu- 
tion of the equations of 10-dimensional IIA supergravity IIA supergrav- 
ity is the dimensional reduction of 11-dimensional supergravity on S^, and 
the D4-brane solution is the dimensional reduction of a M5 brane solution 
that wraps the S^. The near horizon geometry of the M5 brane solution is 
AdSy X S^, and the 11-dimensional equations admit a consistent truncation 



to 7-dimensional equations of motion (1.1) involving gravitational dynamics 



with a negative cosmological constant. Compactifying further on an S^ and 
restricting to the zero momentum sector on this circle, yields a further con- 
sistent truncation of this seven dimensional set of equations. The resulting 
6-dimensional equations are simply the Einstein-dilaton equations about the 
D4-brane background. It follows that the Einstein-dilaton equations consti- 
tute a consistent truncation of the equations of IIA supergravity about the 
D4 near- horizon background (as is easily independently verified). 

It turns out that the fluid dynamics dual to the long wavelength fluc- 
tuations about the thermal M5 brane solution is simply that computed in 
§1.3[ for the special case d = Q. We want to focus on gravitational solutions 
corresponding to fluid flows independent of one of the boundary directions 
x* (the direction of the S^ in the paragraph above). These lie within the 
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six dimensional consistent truncation described above. But these are simply 
the gravitational solutions dual to fluid flows on the world volume theory 
of the D4-brane. In other words, the fluid dynamics of the D4-brane world 
volume theory is a dimensional reduction of the conformal fluid dynamics 
on the world volume of the M5 brane. Moreover the gravitational duals 
to D4-brane fluid flows are very easily obtained from the KK reduction of 



the results of ^.3 Note that the dimensional reduction of conformal fluid 
dynamics results in non-conformal fluid dynamics (e.g. the dimensional re- 
duction of a traceless stress tensor generically has non- vanishing trace). 

It is an interesting and surprising fact that the discussion of the previous 
paragraph generalizes to Dp-branes for all p at the purely algebraic level. 
In every case one can find a consistent truncation of the Einstein-dilaton 
system which, in a purely formal manner, can be regarded as the reduction 
of negative cosmological constant Einstein gravity in a higher (sometimes 
fractional) dimension [28]. This observation immediately yields the fluid de- 
scriptions of arbitrary Dp-brane backgrounds as a dimensional reduction of 



the conformal fluid dynamics derived in ^1.3 



1.6.2 Theories with a deconfinement transition 

So far, we have studied the fluid dynamical description of field theories that 
are 'deconfined' at every temperature, i.e., the free energy is 0{N'^). Con- 
sider, however, a theory like pure Yang Mills at large A^ which undergoes 
a first order deconfinement phase transition at finite temperature. Such a 
system has a dual description in terms of a black hole only above the de- 
confinement temperature; the low temperature phase is given by a gas of 
glueballs, and is thermodynamically indistinguishable from the vacuum at 
leading order in A^ (free energy is 0{1)). 

The Scherk-Schwarz reduction of A/" = 4 Yang Mills on a circle of radius 
R (with anti-periodic boundary conditions for fermions) is a simple exam- 
ple of such a theory. At strong coupling this theory undergoes a first order 
deconfinement transition at T R = 2tt. The gravity dual of the high temper- 
ature phase is simply the S^ compactification of the AdSs planar black hole. 
The gravity dual of the low temperature phase is a so-called AdS-soliton (a 
double analytic continuation of the planar Schwarzschild-AdS black hole, in 
which the role of time and the S^ direction are interchanged). At temper- 
atures much higher than the phase transition, the effective 3-dimensional 
low energy theory is simply the dimensional reduction of the 4-dimensional 



conformal fluid system derived in earlier subsections (just as in ^1.6.1) 



However, at the phase transition temperature we have a new phenomenon; 
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there exists a new static solution of the equations of motion, a co-dimension 
one domain wall that interpolates between the AdS-soliton and the S^ com- 
pactification of the planar AdSs black hole. Unfortunately this solution has 
been constructed only numerically [29j. The domain wall is static in these 
solutions because of a pressure balance on the two sides (recall that the free 
energy, and hence the pressure, of the two phases are equal at a phase transi- 
tion temperature). This configuration is that of a fluid with a boundary; the 
effective low energy fluctuations of this system consist of boundary modes 
(like waves on the surface of water) in addition to the bulk modes discussed 
so far in this chapter. At the ideal fluid level the action for boundary de- 
grees of freedom is captured by a single parameter, the surface tension of 
the boundary (computed from the numerical solution). 

Already using the ideal fluid action including boundary terms, it has 
proved possible to construct many stationary solutions of the fluid equations. 
These solutions, called plasma-balls and plasma-rings, have dual descriptions 
as black holes, black rings, and (in higher dimensions) black objects of more 
exotic topology |30j . The effective action for surface degrees of freedom has 
not been worked out at higher orders in the derivative expansion, and ap- 
pears to be an interesting exercise. 

One particularly interesting static solution of the equations of ideal fluid 
dynamics including boundary terms is the plasma-tube; a configuration con- 
sisting of a domain wall that interpolates from the vacuum to the high tem- 
perature phase at the phase transition temperature, followed by a second 
parallel domain wall at separation L that interpolates back to the vacuum. 
Such a fluid configuration is the two dimensional analogue of a 3-dimensional 
cylindrical tube of fluid, and is well known to undergo a famous fluid dy- 
namical instability (to droplet formation) called the Rayleigh instability. For 
real fluids such as water, the endpoint of the Rayleigh instability is a series 
of disconnected droplets. Now the gravitational dual of the plasma tube is 
an inflnitely long black string in 5-dimensional gravity. This solution has the 
well known Gregory-Laflamme instability [Ch.GL] which, apparently, is dual 
to the Rayleigh instability in the long wavelength limit. The boundary dual 
of a series of disconnected droplets, on the other hand, is a series of discon- 
nected black holes. This discussion at least strongly suggests that the end 
point of the Gregory-Laflamme instability consists of localized black holes 
|31j . Note that the fluid description breaks down near the 'pinch off' point; 
the actual description of topology change in this process requires the use of 
the full fleld theory (e.g. details of interactions between water molecules in 
the case of water) . 
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1.6.3 Charged fluids and anomalies 

Under the AdS/CFT correspondence a global symmetry in the boundary 
maps to a gauge symmetry in the bulk. This suggests that there should 
be a duality between the long wavelength asymptotically AdS planar black 
hole solutions of the Einstein-Maxwell theory (with negative cosmological 
constant) and the equations of charged fluid dynamics. This is a useful ex- 
tension as fluids of interest in experimental situations conserve one or more 
C/(l) charges in addition to energy and momentum. For instance, the flow 
of air in the atmosphere conserves air molecule number. 



It is conceptually straightforward to generalize the set-up of section ^1.3 
to the study of locally thermalized charged planar black holes. The starting 
point is the construction of spacetimes that tubewise approximate Reissner- 
Nordstrom AdS black hole solutions with locally varying temperature, chem- 
ical potential and velocity. A perturbation expansion entirely analogous to 
the one outlined in §1.3| then constructs the gravitational solutions to the 
Einstein-Maxwell-Chern-Simons theory (which forms a consistent trunca- 
tion of IIB supergravity on AdS5 x S^) dual to charged fluid flows order by 
order in the derivative expansion j32^ l33] . This procedure also determines 
the equations of charged fluid dynamics to order by order in the derivative 
expansion. The results of this analysis turned out to throw up a surprise 
purely from the viewpoint of charged first order fluid dynamics, as we now 
explain. 

The equations of charged fluid dynamics are the conservation of the charge 
current 

Var = (1.46) 



together with the conservation of the stress tensor (1.3). Concrete fluid dy- 
namical equations require constitutive relations that express the field redef- 
inition invariant parts of the stress tensor and charge current in terms of 
expressions of first order in the derivative of fluid fields. The charge current 
for charge density q takes the form 

r = qu'' + JSiss, (1-47) 

where J^iss represents the contribution of terms with one or more derivatives 
of the fluid fields to the charge current, and is to be viewed as the charge 
current analogue of 11" ; similarly to (|1.6|) we expand 



71=1 

Standard textbook analyses assert that the most general allowed form at 
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first order for the constitutive relations of a relativistic charged fluid are the 



first equation of (1.7) along with 



^''c (j(i) + J^i^b^li)^) =^K, V,^ ^ -P"V,^ + ^p^ (1.49) 



while the second equation of (1.7) gets modified to 
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In (1.49) F is the non-dynamical background electromagnetic field that 



couples to the U{1) current J" in (1.46) and // is the chemical potential of 
the fluid. Provided that 

r/>0, K>0, /3>0, 

these expressions are consistent with the positivity of divergence of the fluid 
entropy current 
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However it was found by explicit computation that the fluid dual to the 
asymptotically AdS-Einstein-Maxwell-Chern-Simons system has constitu- 



tive relations that differ from those of (1.49) in the following fashion: the 



r.h.s. of (1.49) includes new terms proportional to fluid vorticity oj"" and rest 



frame magnetic field B'^ where 



1 



OJ 



.abed 



UbdcUd , 



B" 



-e-'^'utF.d. 



In a beautiful paper [3l] pointed out the reason for the appearance of these 
new terms. When the U{1) current has a global U{1)^ triangle anomaly (as 
is true of the field theory dual to a bulk system with a 5-dimensional Chern- 



Simons term), (1.52) has an additional term on its r.h.s. proportional to this 
anomaly. This term spoils the positivity of the divergence of the canonical 
entropy current in the presence of such a field. It is however consistent with 
the positivity of the divergence of a modified entropy current provided that 



modifications are also made to the r.h.s. of ( 1.49 ) . More concretely, positivity 



of the entropy current in every conceivable circumstance requires that, in 
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addition to the first equation of (1.7) and (1.50) 



Js = Jean + crij(^"' + CFB B" 



^" {jti) + J^i^bIl\l^))='^Vi+~^^^"+~^BB- 
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(1.54) 



t2 



V 2p„ + P J 2 P + P 

and fci and k2 are integration constants. Furtlier imposition of CPT invari- 
ance forces k2 to vanish. 

This explanation accounts for the additional transport coefficients in the 
AdS/CFT duahty, but apphes more generally to every fluid flow with a 
U{\Y anomaly. The effect of these new transport coefficients may well turn 
out to have experimentally measurable effects in the relativistic heavy ion 
collisions or in neutron or quark stars. 



1.6.4 Holographic superfluid hydrodynamics 

It was pointed out by [35] that charged asymptotically AdSs planar black 
holes are sometimes unstable in the presence of charged scalar fields. The 
endpoint of this instability is a hairy black hole: a black hole immersed in 
a charged scalar condensate. The AdS/CFT correspondence maps the hairy 
black hole to a phase in which a global U{1) symmetry is spontaneously 
broken by the vacuum expectation value of a charged scalar operator (see 
[Ch.CM] for further discussion). In condensed matter physics a phase with 
a spontaneously broken global U{1) symmetry is referred to as a superfluid. 
The variables of relativistic superfluid dynamics consist of two velocity 
fields, the normal fluid velocity u"" and a superfluid velocity field nf , together 
with a temperature and chemical potential field. The superfluid velocity is 
the unit vector in the direction of —^a where ^a is the gradient of the phase of 
the scalar condensate. Conservation of the stress tensor and charge current 
together with the assertion that ^a is curl free constitute the equations of 
superfluid dynamics. These equations form a closed dynamical system once 
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they are supplemented with constitutive relations that express the stress 
tensor, charge current and the component of S^a along the normal velocity, 
as functions of the fluid dynamical variables. 

It has proved possible to apply the fluid/gravity map to hairy black holes 
to derive the constitutive relations for holographic superfluids, with interest- 
ing results. The theory of perfect superfluids was worked out by Landau and 
Tisza in the 1940s. In a beautiful recent work |36] have used the equations 
of Einstein gravity to re-derive Landau-Tisza equations for superfluids that 
admit a holographic description. The theory of first order dissipative correc- 
tions to the equations of Landau-Tisza superfluidity was most completely 
spelled out in [37J. Calculations done within the fluid/gravity framework 
have led to the realization that the 13 parameter Clark-Putterman equa- 
tions derived therein miss one parameter (under the assumption of parity 
invariance for the superfluids) and 6 more parameters (if the superfluids are 
not assumed to preserve parity). A completely satisfactory framework for 
superfluid hydrodynamics has been developed only very recently [7], and 
the fluid/gravity map has played a major role in this development. 



1.7 Relation to other developments 

Having surveyed the fluid/gravity correspondence and its various applica- 
tions, we finally describe connections with other approaches. 

1.7.1 Implications for Israel- Stewart formalism 

One useful application of the fluid/gravity correspondence is an improve- 
ment on the 'causal relativistic hydrodynamics', also known as the Israel- 
Muller- Stewart formalism j38l I39| . To appreciate the context, recall that 
a conventional theory of relativistic dissipative (i.e., irreversible) hydrody- 
namics, which is first order in time derivatives, is described in terms of a 
parabolic system of differential equations, leading to instantaneous propa- 
gation of signals. While these apparently a-causal modes lie outside of the 
long wavelength regime of validity of the hydrodynamical formulation as 
discussed in [lO], they nevertheless lead to conceptual and computational 
problems. To capture the dissipative physics, [38] observed that second order 
terms are needed in the entropy current. These render the system hyperbolic, 
thereby providing a good initial value formulation. However, the particular 
terms added are not all the possible ones consistent with the symmetries, so 
as such, the construction is somewhat ad-hoc. Indeed, [T6] observed in the 
context of conformal fluid, that the terms added do not maintain conformal 
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invariance of the system, manifesting the incompleteness of the approach. 
The fluid/gravity construction in effect prescribes the correct completion 
to render the full system causal, as well as manifestly consistent with the 
symmetries. We expect that due to the gravitational dual, causality will be 
guaranteed at all orders in the derivative expansion. 



1.7.2 The black hole membrane paradigm 

Perhaps the most salient feature of the fluid/gravity correspondence is the 
fact that the horizon dynamics (which in this case prescribes the dynamics 
of the entire spacetime) is governed by hydrodynamics. At the face of it, 
such type of relation is not new; in fact for several decades relativists have 
explored the idea that spacetime, or important aspects thereof like black 
hole horizons, might resemble a fluid. Early indications include black hole 
thermodynamics [121 [H] developed in the 70's, analog models of black holes 
|41j initiated in the early 80's, and most strikingly the black hole Membrane 
Paradigm [42^143] formulated in the late-70's. The latter realizes the idea that 
for external observers, black holes behave much like a fluid membrane, en- 
dowed with physical properties such as viscosity, conductivity, and so forth. 
In particular, the dynamics of this membrane is governed by the familiar 
laws of fluid dynamics, namely the compressible Navier-Stokes equations. 

Motivated by the superficial similarity between the membrane paradigm 
and the fluid/gravity correspondence, recently [Ml SS] have attempted to 
formulate a precise derivation of the former. In [H] Einstein's equations 
in the bulk are projected onto a null hypersurface and then expanded in 
gradients along the hypersurface. On the other hand |45j show that one 
can systematically find a solution to vacuum Einstein's equations which 
describes the near-horizon geometry of a generic non-degenerate black hole 
in the long wavelength regime. 

Within the fiuid/gravity correspondence, the entire spacetime evolution 
is mapped to the dynamics of a conformal fluid, which, albeit reminiscent of 
the membrane paradigm, has one important twist: the membrane lives on 
the boundary of the spacetime (which is unambiguously defined and admits 
a fluid description with well-defined dynamics), and gives a perfect mirror 
of the full bulk physics. This "membrane at the end of the universe" pic- 
ture is a natural consequence of the holographic nature of the fluid/gravity 
correspondence. 
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1.7.3 Blackfolds 

As in the fluid/gravity correspondence and the membrane paradigm type 
ideas, the blackfold approach to constructing higher-dimensional black holes 
(discussed in [ch.BF] of this book) likewise asserts that the effective theory 
describing the long wavelength dynamics of black hole horizons can be ex- 
pressed in terms of fluid dynamics. However, there are several important 
differences between these descriptions. Since the blackfold 'fluid' pertains 
to the effective world-volume dynamics of an extended black object as seen 
from far away, the intrinsic dynamics typically has to be supplemented by 
extrinsic dynamics, describing how the blackfold embeds in the ambient 
spacetime. On the other hand, in the fluid/gravity correspondence the fluid 
resides on the boundary of asymptotically (locally) AdS spacetime, so there 
is no issue with extrinsic dynamics. 

Moreover, although the blackfold formalism is most naturally formulated 
in asymptotically flat spacetime, by a suitable separation of scales, one 
can in principle consider blackfolds with any asymptotics. In contrast, the 
fluid/gravity correspondence concerns asymptotically AdS black holes. On 
the other hand, unlike all the above-mentioned approaches, fluid/gravity 
is the only one where there is a known physical microscopic origin to the 
fluid: it is the effective behavior of the dual field theory residing on the AdS 
boundary. 



1.8 Summary 

The fluid/gravity correspondence provides a natural way to map solutions 
of fluid dynamics into those of gravity, enabling one to construct time- 
dependent, inhomogeneous black hole solutions to Einstein's equations, re- 
taining full non-linearity. An interesting aspect of the construction is the 
manner in which classical gravity can be moulded to fit naturally with ef- 
fective field theory intuition to extract approximate solutions. While the 
construction itself arose from the gauge/gravity correspondence, it is clear 
that it can be implemented in greater generality. 

Apart from providing interesting insights into the dynamics of gravity, 
the map has played an important role in clarifying various issues in fluid dy- 
namics. The role of quantum anomalies in hydrodynamical transport, and 
generalizations of fluid dynamics to systems with spontaneously broken sym- 
metries, are two examples where the fluid/gravity map has served to eluci- 
date the underlying physics cleanly. The physical points seem much simpler 
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to understand from the gravitational perspective; aided by this intuition one 
can re-evaluate the hypotheses of traditional descriptions of fluids. 

The fluid/gravity map suggests several extremely interesting technical, as 
well as conceptual, questions for the future. Some of these are: Does the 
gravitational viewpoint shed any light on turbulent fluid flows, or questions 
about singularities that develop in finite time from smooth initial data in 
fluid dynamics? Is there a path integral formulation of fluid dynamics at 
finite N, and how does it map to the path integral of bulk gravity? Are the 
corrections to the classical equations of gravity constrained by the require- 
ment of positivity of divergence of an 'entropy current' on an event horizon 
(analogously to, and perhaps even dual to, fluid dynamics)? It seems likely 
that many interesting results remain to be discovered in this general area. 



1.9 Epilogue: Einstein and Boltzmann 

As we have emphasized throughout, the equations of fluid dynamics, for 
which we have an independent field theory intuition, are dual to a long 



wavelength limit of Einstein's equations ( 1.1 ). It is then natural to ask what 



is the field theoretical interpretation of the full dynamical system of equa- 



tions (1.1)? We believe that these equations may be conceptually thought 
of as the strong coupling analogue of (a decoupled sector of) the Boltzmann 
transport equations. 

It is well known that the linearization of the Boltzmann transport equa- 
tions about equilibrium yields an infinite set of 'quasinormal modes', i.e., 
solutions to the equations of motion that all decay to zero (returning the 
system back to equilibrium) at late times. Exactly d of these quasinor- 
mal modes are massless (in the sense that they are static in the infinite 
wavelength limit). In textbooks on statistical mechanics, fluid dynamics is 
sometimes derived as the non-linear theory of this finite set of Boltzmann 
'quasinormal modes'. The remaining quasinormal modes are 'fast modes' 
that decay away on a time scale set by the mean free path of kinetic theory. 

Similarly, the linearization of the equations of gravity about the planar 
black hole has d massless quasinormal modes and an infinite set of mas- 
sive quasinormal modes. In direct analogy with the work on the Boltzmann 
equations, the fluid/gravity correspondence constructs the equations of fluid 
dynamics as the non-linear theory of these massless modes (effectively by 
integrating out the massive modes, order by order). For this reason it is 
natural to think of the full set of Einstein's equations in the presence of an 
event horizon (including all quasinormal mode degrees of freedom) as the 
strong coupling analogue of the Boltzmann transport equations. 
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An important property of the Boltzmann transport equations is that they 
are irreversible; they obey the Boltzmann i^-theorem (which asserts that a 
certain functional of kinetic variables called H always increases in time and 
is maximum in equilibrium). In direct analogy, Einstein's equations, together 
with the assumption of regularity of future event horizons (and physical en- 
ergy conditions), always obey the classic area increase theorem of general 
relativity. This suggests that the better analogy is between the Boltzmann 
transport equations and Einstein's equations plus the condition of regular- 
ity of the future event horizon. The last condition breaks the time reversal 
invariance of Einstein's equations. In fact, the requirement that the future 
event horizon stay regular was a crucial element in our implementation of 
the fluid/gravity map. The Boltzmann theorem has a local analogue in fluid 
dynamics; it maps to the statement that the equations of fluid dynamics are 
accompanied by a local entropy current that whose divergence is everywhere 
non-negative. The area increase theorem of general relativity can be used 
to construct such an entropy current for the fluid dynamics generated from 
the fluid/gravity map. 

Just like the Boltzmann equations, the system of gravitational equations 



(1.1) can be used to study the approach to equilibrium from a highly non 
equilibrated starting point. In simple studies of equilibration using Einstein's 
equations |461 I47j . the equilibration time, measured by the time taken for 
fluid dynamics to take over as the effective description, turns out to be 
extremely rapid. In more complicated situations the equilibration process 
displays sharp phase transitions associated with Choptuik phenomena. In- 



deed the equations (1.1) undoubtedly contain a host of dynamical delights 
for the intrepid gravitational and statistical physicist; it seems clear that 
the fluid/gravity map is merely the tip of an iceberg's worth of connections 
between gravity and statistical physics. 
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